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Resumen 


En este trabajo se analiza el efecto de la resolución espacial de las imágenes 
utilizadas para conformar el espacio Normalized Difference Vegetation Index (NDVD- 
temperatura de superficie (T,). Este espacio es muy utilizado para el monitoreo 
de variables ambientales, tales como la evapotranspiración (ET) y la fracción 
evaporable (FE). El método propuesto por Jiang e Islam (2001) fue seleccionado 
para estudiar el efecto de la resolución espacial en el cálculo de la FE. Para esto, 
se analizaron resoluciones espaciales que varían entre 1 x 1 km, comunes en las 
imágenes de sensores térmicos, hasta 10 x 10 km, consistente con las imágenes 
de sensores de microondas pasivos (MP). Los resultados sugieren que el espacio 
triangular NDVI-T, permite obtener valores razonablemente precisos de FE con 
imágenes de baja resolución espacial, como la de los MP, siempre que el borde 
frío pueda ser determinado con imágenes de resolución espacial 1 x 1 km. Si no es 
posible calcular dicho borde con imágenes de resolución 1 x 1 km, los resultados 


serían inciertos y por tanto deben ser cuidadosamente analizados. 
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Introducción 


El manejo de las actividades agrícola y 
forestal, así como el de los recursos hídricos, 
requiere del balance de energía radiactiva a 
nivel de la superficie terrestre a diferentes 
resoluciones espaciales y temporales. Así, en 
los últimos años se han desarrollado diferentes 
métodos para calcular el flujo de calor latente 
o evapotranspiración (ET) bajo diferentes 
condiciones de humedad de la superficie y 
escalas espacio-temporales (Jackson et al., 1977; 
Jiang e Islam, 2001; Norman et al., 2003; Rivas y 
Caselles, 2004). Estas metodologías requieren de 
datos de temperatura de la superficie terrestre 
(T). Con la tecnología de satélite es posible 
obtener T, con una distribución espacio-tempo- 
ral, que es fundamental para el monitoreo de 
variables ambientales y ecosistemas en general. 
Un ejemplo son los mapas de T. generados 


con información del sensor MODIS (Moderate 
Resolution Imaging Spectroradiometer) y dis- 
tribuidos en forma gratuita a la comunidad 
científica. MODIS es uno de los sensores a bordo 
de los satélites EOS-Terra y EOS-Aqua operados 
por la NASA  (http://modis.gsfc.nasa.gov). 
Estos mapas poseen una resolución temporal 
diaria y una resolución espacial de 1 x 1 km, y 
han sido validados mediante datos de campo 
bajo diferentes condiciones de la superficie (Wan 
et al., 2002; Wan et al., 2004; Coll et al., 2005; Wan, 
2008). 

En los últimos años se han desarrollado 
diferentes metodologías para calcular el 
balance de energía, más precisamente la 
ET, con datos de satélite. Muchos de estos 
métodos se basan en la relación contextual 
entre T, y un índice de vegetación (Jiang e 
Islam, 1999 y 2001; Nishida et al., 2003; Stisen 
et al., 2008). De todos estos métodos, se destaca 
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aquí la simplicidad del método de Jiang 
e Islam (2001), que calcula ET y la fracción 
evaporable (FE) espacialmente distribuida, 
basándose en la interpretación contextual de 
la relación entre la T' y el índice de vegetación 
diferencial normalizado (NDVI) (Jiang e 
Islam, 1999). Básicamente, este método utiliza 
el diagrama NDVI-T, de forma tal que FE 
puede ser calculada con un mínimo número 
de parámetros de calibración. 

La principal ventaja de una formulación 
contextual es que es válida para una gran gama 
de condiciones superficiales y posiblemente 
también sobre un amplio rango de resolucio- 
nes espaciales. Otra ventaja sobresaliente de 
este contexto es que no se requiere de gran 
exactitud en las estimaciones de la tempera- 
tura de la superficie. Jiang e Islam (2001), Jiang 
et al. (2004) y Venturini et al. (2004) demostra- 
ron que la temperatura obtenida con una sola 
banda térmica es suficiente para extraer valores 
razonablemente precisos de FE. Por otra parte, 
Kustas et al. (2003) y French et al. (2005) han 
estudiado la variabilidad de la geometría del 
espacio NDVI-T. a escalas espaciales inferiores 
a 1 km?. Los autores desarrollaron un método 
para bajar la resolución de la T' 

Los sensores térmicos presentan el incon- 
veniente de no poder registrar las ondas 
infrarrojas térmicas bajo condiciones climáticas 
adversas, esto significa que no es posible obte- 
ner valores de Ten días nublados. Una solución 
para este problema serían los sensores de 
microondas pasivos que registran información 
diariamente y operan con longitudes de 
ondas más largas, que no son afectadas por 
las condiciones climáticas ni por la fuente de 
iluminación. Esta fuente alternativa de datos 
podría sustituir los mapas de T, obtenidos con 
sensores térmicos, previo análisis del efecto que 
tiene la mala resolución espacial de los MP. 

En este trabajo se investiga el efecto que 
tiene la resolución espacial de la variable T. en 
el espacio NDVIT, ya que este espacio define el 
cálculo del parámetro 4 del modelo propuesto 
por Jiang e Islam (2001) para la estimación de 
la FE. 


Método de Jiang e Islam 


La metodología desarrollada por Jiang e Islam 
(2001) propone una modificación de la ecuación 
de Priestley y Taylor (1972), reemplazando el 
parámetro Q. de Priestley-Taylor por un nuevo 
parámetro q, como se observa en la ecuación 


(1) 


er=0| 2 [n= 0) 0) 


donde y es la constante psicrométrica; A, la 
pendiente de la curva de presión de saturación 
del vapor de agua; R,, la radiación neta que 
llega a la superficie; G, el flujo de calor del 
suelo, y q, el parámetro de Jiang e Islam. 

La ecuación modificada (ecuación (1)) puede 
ser aplicada a regiones con cobertura vegetal 
heterogénea y con diferentes condiciones de 
humedad en la superficie, a diferencia de la 
ecuación original de Priestley y Taylor, que 
fue desarrollada para superficies saturadas 
(Jiang e Islam, 1999). Esta metodología basa la 
modificación introducida en la interpretación 
del espacio NDVI-T, donde los píxeles adoptan 
una distribución trapezoidal al ser graficados, 
tal como se puede observar en la figura 1 (Jiang 
e Islam, 1999). 
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Figura 1. Espacio contextual NDVI- T, y definición de los 
parámetros del método de Jiang e Islam 
(fuente: Jiang e Islam, 1999). 
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La interpretación física de la relación NDVI- 
T, ha sido objeto de diversas investigaciones y 
discusiones. Estas discusiones se centran en la 
definición del lado superior (o borde cálido) 
del triángulo, en donde los píxeles con altos 
valores de T. tienden a alinearse en forma de 
una recta bien definida en el rango del NDVI 
(ver figura 1). Jiang e Islam (2001) interpretaron 
que el borde superior constituye el mínimo 
valor de ET para cada clase de NDVI, mientras 
que el borde frío (o lado inferior del triángulo) 
representa la máxima ET. De esta forma, la 
magnitud de ET para cada píxel varía dentro de 
los límites del triángulo (ver figura 1). Por esto 
es necesario obtener para cada píxel el valor 
de q, el cual se calcula mediante una interpola- 
ción lineal entre ambos lados del triángulo. 

El parámetro $ puede variar desde cero 
para un suelo desnudo y seco hasta 1.26, para 
una superficie saturada. Cabe aclarar que 
el término superficie saturada se considera 
de ahora en más como el conjunto suelo- 
vegetación saturado de agua o un píxel donde 
hay una lámina de agua en superficie. La 
variación de q depende de las características 
físicas de la superficie. Estas características 
estarían asociadas con humedad, conductancia 
y temperatura de la superficie (Jiang e Islam, 
2001). 

Para obtener el valor de ( se debe definir la 
pendiente del lado superior del triángulo para 
los píxeles cuyos NDVI estén comprendidos en 
el rango [0, 1]. El valor mínimo global de 4 (4,.... 
= 0) está asociado con una superficie seca sin 
vegetación, con un NDVI=0 y una temperatura 
alta denominada T'., que se obtiene como la 
intersección del borde superior del triángulo 
con la recta NDVI = 0. El valor máximo 
global de 6 (6. = 0. = 1.26) se define para la 
temperatura media de los píxeles en los que 
se observe agua libre (T, .). Una vez definidas 
Tu Y Lay Correspondientes a 4 =0 y 4 = 1.26, 
respectivamente, se realiza la interpolación 
lineal de la que resulta la fórmula siguiente 
(Venturini et al., 2004): 


TT má 
En a 1 máx 2 
i La e El : 


donde T, es la temperatura del píxel i y 4, = 
1.26. 


Tratamiento de las imágenes 


La temperatura de cualquier superficie posee 
una variabilidad temporal notable, con dife- 
rencias que pueden llegar a los 10 “K entre 
el día y la noche. La variabilidad espacial de 
T. depende del tipo de cobertura vegetal que 
caracteriza el área en estudio. Se considera que 
la temperatura de un píxel es el resultado de la 
cantidad de radiación emitida por una superfi- 
cie mixta (suelo y vegetación), la que es captada 
por el sensor como una señal única (Venturini 
et al., 2008). Teniendo en cuenta esto, se centró 
la atención en las imágenes de T, debido a la 
dinámica que dicha variable presenta frente al 
NDVL, el que puede considerarse casi constante 
durante 10 o 15 días (Huete ef al., 1999; Huete 
et al., 2002). 

Con el fin de obtener los espacios NDVIT, 
se procesaron imágenes de T, correspondien- 
tes al producto MOD11 (http: / /modis-land. 
gsfc.nasa.gov /temp.htm). Se georregistraron 
las imágenes a partir de las coordenadas 
geográficas (latitud y longitud) originales de las 
mismas. Posteriormente, se procedió a cortar 
el área que comprende la región en estudio. 
Luego, las imágenes de T fueron agregadas 
mediante simple promedio aritmético a partir 
de la imagen original a resolución 1 x 1 km. 
Por ejemplo, para obtener T, a resolución 2 x 
2 km, se promediaron de a cuatro los píxeles 
a partir de la imagen de T, a resolución 1 x 1 
km. Seguidamente se calculó el NDVI a partir 
de las bandas del rojo e infrarrojo cercano. 
Las imágenes de NDVI se trabajaron con su 
resolución original (1 x 1 km), a fin de conservar 
el detalle. Como el NDVI es obtenido en fun- 
ción de dos bandas, una en el espectro visible 
y la otra en el infrarrojo cercano, el proceso de 
agregación de este índice requiere de un análi- 
sis que está fuera del alcance de este trabajo. 
Por lo tanto, para poder equiparar el tamaño de 
las matrices de NDVI con las de T' se recurrió al 
“pixelado” de estas últimas. Se subdividieron 
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los píxeles agregados en unidades menores, 
más específicamente en píxeles de 1 x 1 km, 
de igual valor. Por ejemplo, cada píxel de 
resolución 2 x 2 km fue dividido en cuatro 
píxeles de 1 x 1 km de igual temperatura. Este 
procedimiento se denominó “pixelado” de las 
imágenes de T, agregadas. 

Con las matrices de T, y NDVI de igual 
dimensión se obtuvieron los espacios contex- 
tuales NDVI-T. y se compararon las geometrías 
de los diagramas obtenidos. Finalmente, se 
obtuvieron las imágenes de ( y se contrastaron 
los resultados para las diferentes resoluciones 
analizadas. De esta forma se simuló la fusión 
de imágenes de T| provenientes de un sensor 


MP con imágenes de NDVI, sin degradar la 
resolución espacial de estas últimas. 


Área de estudio 


El área de estudio corresponde a la cuenca 
inferior del río Salado, provincia de Santa 
Fe, Argentina, ubicada entre los 29 y 31” de 
latitud sur, y los 62” y 60.5” de longitud oeste 
(figura 2). 

El río Salado se extiende a lo largo de dos 
mil kilómetros. Nace en la provincia de Salta 
y desemboca en el sistema del río Paraná. La 
cuenca hidrográfica comprende unos 58 000 
km?, de los cuales unos 30 000 km? aportan 
directamente al río Salado. 


Mapa de NDVI 


529” W62" 


S31* W60.5* 


Figura 2. Ubicación geográfica del área de estudio, Santa Fe, Argentina. 


Girolimetto ef al., Influencia de la resolución espacial de los diagramas NDVI-T en el cálculo de la fracción evaporable 


En los últimos treinta años, la cuenca 
experimentó cambios hidrológicos, como el 
aumento del caudal medio del río debido al 
incremento de precipitaciones (FICH-INA- 
INTA, 2007). Paralelamente a los cambios 
de las variables externas al sistema, se han 
introducido cambios en el uso del suelo, 
creciendo la superficie cultivada en detrimento 
de la utilizada para ganadería. La agricultura 
es intensiva en este área, mas precisamente el 
cultivo de oleaginosas (soja, girasol y maíz) es 
el predominante. Paulatinamente, la dinámica 
hídrica natural también fue modificándose 
por la construcción de obras de canalización, 
alcantarillado y puentes. Laimportancia de esta 
región radica en el hecho de que la provincia de 
Santa Fe es una de las principales productoras 
de soja (http: / /www.santafe.gov.ar/). 


Datos utilizados 


A fin de disminuir la incertidumbre que gene- 
raría el uso de mapas no validados derivados 
de MP, se simularon las imágenes de Ta 
diferentes resoluciones espaciales. Para esto se 
utilizaron imágenes obtenidas con un sensor 
térmico. La metodología utilizada para obtener 
las imágenes de T, a resolución de 1 x 1 km ha 
sido ampliamente validada (Wan y Dozier, 
1996). Esto significa que los errores de las 
imágenes de T están acotados y son conocidos, 
no así los derivados de sensores MP. 

El programa EOS (Earth Observing System), 
desarrollado por la NASA, consiste en una 
serie de satélites que orbitan la Tierra en 
forma sincronizada. El primer satélite del 
programa EOS, EOS-Terra, fue puesto en órbita 


en diciembre de 1999, con el sensor MODIS 
(Moderate Resolution Imaging Spectroradiometer) 
entre los cinco sensores a bordo. EOS-Terra 
orbita a una altitud de 705 km y recorre la 
Tierra varias veces al día. La órbita del satélite 
va de norte a sur y cruza el ecuador por la 
mañana. 

MODIS es uno de los instrumentos a 
bordo de EOS-Terra. MODIS proporciona 
información en 36 bandas espectrales en el 
rango de longitudes de onda de 0.4 — 14.385 
um (visible a infrarrojo térmico), con una 
resolución espacial que varía de los 250 a los 
1 000 m, y un ancho de barrido de 2 330 km. 
En la actualidad se distribuyen 44 productos 
derivados de las 36 bandas, los que se dividen 
en las siguientes categorías: Calibración, 
Atmosféricos, Terrestres, Criósfera y Oceánicos 
(http: / /modis.gsfc.nasa.gov). 

En este trabajo se utilizaron los productos 
MOD02 y MODl1. El producto MODO2 
contiene la radiancia geo-registrada de las 36 
bandas. Las imágenes correspondientes a la 
banda 1 (rojo) y la banda 2 (infrarrojo cercano) 
fueron utilizadas para obtener el índice de 
vegetación NDVI. El producto MOD11 ofrece 
la temperatura de superficie georreferenciada 
y los valores de las emisiones de las bandas 
térmicas 31 y 32 para píxeles de 1x1y5x5 
km, respectivamente (Wan y Dozier, 1996). 
Las imágenes de T, del producto MOD11 se 
utilizaron para trazar el espacio NDVIT, y 
para obtener las imágenes del parámetro 
(Venturini et al., 2004; Batra et al., 2006). 

Se procesaron tres imágenes MOD11 y 
dos MODO?, correspondientes al mes de sep- 
tiembre de 2004. La fecha y hora de adquisición 


Cuadro 1. Año de adquisición de las imágenes, día del año, producto MODIS, hora de adquisición 
y calidad de las imágenes. 


mn E e Hora de adquisición Calidad de la imagen 
Año Día del año Producto 
(UTC) (% nubes) 
2004 259 MODO02 MOD11 14:15 0-2 
2004 266 MOD11 14:20 0.8 
2004 268 MODO02 MOD11 14:10 0-9 


Ciencias del Agua, vol. Il, núm. 3, julio-septiembre de 2011 


Tecnología y 


Ciencias del Agua, vol. IL, núm. 3, julio-septiembre de 2011 


« vj30]0u32], 


Girolimetto ef al., Influencia de la resolución espacial de los diagramas NDVET en el cálculo de la fracción evaporable 


de cada una de las imágenes y el porcentaje de 
nubes observado se presentan en el cuadro 1. 


Resultados 


Se agregaron los píxeles de las imágenes de 
MOD11 por promedio aritmético, obteniéndo- 
se imágenes con resoluciones espaciales de 
2x2,4x4,8x8 y 10 x 10 km. Con el fin de 
analizar los cambios regionales de T, a medida 
que se degrada la resolución espacial, se 
calcularon los valores de T, mínimos, medios 
y máximos de las imágenes agregadas. El 
análisis de los histogramas de frecuencia de 
las imágenes agregadas evidenciaron que para 
todas las resoluciones estudiadas, las máximas 
y mínimas frecuencias, aproximadamente, se 
registran para los mismos valores de T.. Estos 
resultados no se muestran en este trabajo, pero 
pueden ser consultados en Girolimetto et al. 
(2008). 

Si bien la información térmica media 
regional no cambiaría con la disminución de 
la resolución espacial, esto no implica que la 
incertidumbre de las estimaciones hidrológicas 
derivadas de diferentes imágenes agregadas 
se mantenga (Venturini et al., 2004; Batra et al., 
2006). 

La calidad de las imágenes de T', analizada 
en términos de tamaño de píxel, afecta direc- 
tamente la construcción del espacio triangular 
NDVLIT. y, por lo tanto, impacta sobre todas 
las variables ambientales que se deriven del 
mismo. Por ejemplo, el parámetro 4 de Jiang e 
Islam se obtiene con los espacios NDVI-T.. En 
la figura 3 se observan las variaciones que se 
producen en la nube de puntos al aumentar el 
tamaño del píxel, evidenciándose la pérdida de 
los valores extremos a medida que disminuye 
la resolución espacial. 

En el caso particular del método de 
Jiang e Islam, los bordes inferior y superior 
del triángulo son los elementos clave para 
la determinación de q. Es por ello que se 
analizaron por separado las particularidades 
de cada uno de dichos bordes. 


Análisis del borde superior del espacio NDVI- 
Í, 

Las características del borde superior de 
los diferentes espacios triangulares están 
resumidas en el cuadro 2. Se observan 
diferencias en las pendientes de los límites 
superiores de los triángulos NDVI-T.. Dichas 
diferencias se traducen en variaciones del 
parámetro T',. del modelo de Jiang e Islam 
(para NDVI = 0) entre los 2 y 4 “K. Se destaca 
que en la búsqueda bibliográfica no se han 
encontrado estudios sobre la variabilidad 
del parámetro T,.. Sin embargo, Kustas et 
al. (2003) y Venturini et al. (2004) encontraron 
diferencias similares al analizar la temperatu- 
ra radiométrica de una banda térmica para 
diferentes sensores y a resoluciones espaciales 
de hasta 1 x 1 km. En general, estos resultados 
indicarían que imágenes de T, con píxeles 
más grandes de 1 x 1 km no influirían en el 
diagrama triangular, más específicamente en 
el valor de T',. correspondiente a NDVI=0. 


Análisis del parámetro q, manteniendo el 
borde inferior T_. constante 


Se entiende por T”, . constante a la temperatura 
media de aquellos píxeles con NDVI < 0 
obtenida para la resolución espacial original, en 
este caso de 1 x 1 km. El T,, así determinado se 
mantiene constante para las demás resolucio- 
nes espaciales y representa el borde inferior 
del triángulo para todos los espacios NDVI-T, 
obtenidos. 

Las imágenes de q que se muestran en 
la figura 4 se generaron con la ecuación (2), 
con T,, expuestas en el cuadro 2 y con T,... 
constante. Los valores medios y desvíos 
estándar de ( se resumen en el cuadro 3. Como 
puede observarse, no se identifica un patrón 
de variación de $ a medida que se cambia 
la resolución espacial de T. Las diferencias 
encontradas en el valor medio regional de q 
son inferiores a 0.09 (14% del valor medio). Los 
resultados son compatible con lo observado por 
Batra et al. (2006). Estos autores determinaron 
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Figura 3. Ejemplo de los diagramas NDVI-T.. Día 268 /2004 para resoluciones: (a) 1 x 1 km, (b) 2 x 2 km, (c) 4x 4 km, 
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(d) 8x 8 km, (e) 10 x 10 km. 
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Cuadro 2. Características de los espacios triangulares para 
diferentes resoluciones espaciales. 


2 Da A año 00 
Resolución NOV a, Pendiente Trá 10 
1 km 0.68 -0.0518 310.83 
2 km 0.68 -0.0667 309.91 
4 km 0.68 -0.0603 311.26 
8 km 0.68 -0.0751 309.62 
10 km 0.68 -0.0608 311.06 
Día 266, año 2004 
Resolución NDVI Pendiente Tas 10 
1 km 0.66 -0.0699 312.16 
2 km 0.66 -0.0506 315.00 
4 km 0.66 -0.0573 313.88 
8 km 0.66 -0.0701 S1237 
10 km 0.66 -0.0633 313.05 
Día 268, año 2004 
Resolución NOV a, Pendiente Tirá, 10 
2 km 0.66 -0.0329 327.50 
4 km 0.66 -0.0364 SIDO 
8 km 0.66 -0.0489 322.37 
10 km 0.66 -0.0470 321.85 


el valor de 4 para diferentes sensores con 
resolución espacial 1 x 1 km, obteniendo 
diferencias de alrededor del 10% en el valor 
medio regional de q. Venturini et al. (2004) 
encontraron diferencias en el valor medio de 
4 de aproximadamente un 15% para diferentes 
sensores, siempre con una resolución 1 x 1 km. 

El desvío estándar (S) de b es inferior al 20% 
de la media, lo que estaría indicando una baja 
dispersión de ( en torno al valor medio. En 
otras palabras, la región presenta valores de FE 
relativamente homogéneos. 


Los resultados estarían demostrando que 
el método desarrollado por Jiang e Islam daría 
resultados aceptables en términos de ( para el 
área en estudio, con imágenes cuyos píxeles 
varían de 1 x 1 km a 10 x 10 km, siempre que 
T,,, $e calcule con la imagen de mejor resolu- 
ción espacial. 


Análisis del parámetro q con el borde inferior 
T,,,, Variable 

Se entiende por T'. variable a la temperatura 
media de aquellos píxeles con NDVI < 0, 
obtenida para cada una de las imágenes 
agregadas, la que representa el borde inferior 
del triángulo NDVI-T. Cabe destacar que el 
problema de píxeles mixtos (píxeles donde 
podrían coexistir diferentes tipos y usos de 
suelos y agua libre) es mayor a medida que el 
tamaño del píxeles más grande. Esto implicaría 
que difícilmente se encuentren píxeles puros, 
con agua en superficie. 

En este análisis se trazaron los triángulos 
con los bordes cálidos analizados en el apartado 
“Análisis del borde superior del espacio NDVI- 
T/”; pero en este caso, T', . se calculó para cada 
imagen agregada. Se observó que los triángulos 
se abaten a medida que la resolución espacial 
disminuye. Esto fue consistente en todos 
los días analizados. El cálculo de ú se realizó 
nuevamente con la ecuación (2), con los valo- 
res de T,.. del cuadro 2, pero ahora con T,.. 
estimadas para cada resolución espacial. La 
figura 4 muestra las imágenes de q generadas a 
partir de estos triángulos (ver espacios NDVI-T. 
en la figura 3). 

En el cuadro 4 se muestran los valores 
medios y desvíos estándar de las imágenes 


Cuadro 3. Comparación de las medias y desvíos estándar (S) de los valores de ( para diferentes resoluciones 


para T,,, constante. 
mín 


z 1x1km 2x2km 4x4 km 8x8km 10 x 10 km 
E Media S Media S Media S Media S Media S 
259/04 0.62 0.061 0.58 0.069 0.61 0.063 0.56 0.075 0.60 0.068 
266/04 0.49 0.092 0.56 0.079 0.53 0.084 0.49 0.098 0.51 0.099 
268/04 0.63 0.083 0.67 0.074 0.64 0.083 0.58 0.096 0.58 0.089 
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Mapa de q Mapa de q 
MOD11 -4x 4 km MOD1 - 4 x4 km 
268/2004 268/2004 


A A 
0.38 0.51 0.67 0.80 N/D 0.38 0.51 0.69 0.80 N/D 
a) c) 
Mapa de q Mapa de q 
MOD11 -4x4km MOD11 - 4x4 km 
268/2004 268/2004 


0.20 0.36 0.52 0.71 N/D 0.22 0.46 0.70 0.94 N/D 
d) e) 


Figura 4. Imágenes de 4 — Día 268/2004 para resolución 4 x 4 km (a) y 10 x 10 km, (b) con T,, constante y 4 x 4 km, 
(c) 10 x 10 km (d), para T,, variable. Nótese la diferencia en las escalas numéricas. 


Cuadro 4. Comparación de las medias y desvíos estándar (S) de los valores de $ para diferentes 


resoluciones para T_.. variable. 
mín 


. 1x1km 2x2km 4x4 km 8x8km 10 x 10 km 
me Media S Media S Media S Media S Media S 
259/04 0.62 0.061 0.59 0.071 0.67 0.069 0.65 0.087 0.67 0.075 
266/04 0.49 0.092 0.62 0.087 0.63 0.099 0.78 0.136 0.76 0.130 
268/04 0.63 0.083 0.72 0.081 0.89 0.106 0.97 0.129 0.96 0.139 
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de ( para los diferentes días y resoluciones 
espaciales analizadas. En general, el valor 
medio de Q para cada uno de los días se 
incrementa a medida que aumenta el tamaño 
del píxel de T, debido a que se va perdiendo 
la definición de los valores extremos de T.. Se 
observa que las diferencias en el valor medio 
estimado de 4 oscilan entre 0.03 y 0.34 (de 
5 a 54% del valor medio), evidenciándose 
las mayores a partir de tamaños de píxeles 
superiores a los 4 x 4 km. 

Los resultados obtenidos indicarían que el 
método desarrollado por Jiang e Islam arroja- 
ría resultados de q aceptables con imágenes, 
cuya resoluciones espaciales varíen entre 1 
x 1 y 4 x 4 km. Para resoluciones espaciales 
superiores a los 4 x 4 km, los resultados del 
método diferirían de los obtenidos con las 
imágenes a resolución original, por lo que se 
debe tener precaución al analizar los resultados 
en estos casos (Venturini et al., 2004; Batra et al., 
2006; Girolimetto et al., 2007). 


Conclusiones 


En la actualidad, los sensores remotos ofrecen 
información en un amplio rango del espectro 
electromagnético que permite estudiar pro- 
piedades de la superficie y ecosistemas en 
general. La información proveniente de las 
diversas misiones satelitales ha sido utilizada 
por la comunidad científica para el desarrollo 
de diferentes metodologías, como las publi- 
cadas por Norman et al. (1995), Jiang e Islam 
(2001), Nishida et al. (2003), Rivas y Caselles 
(2004), y Venturini et al. (2008), entre otros. La 
interpretación física del espacio NDVET, ha 
sido objeto de diversos estudios. La mayoría 
de ellos coinciden en que el espacio trapezoidal 
resume el balance energético de la región bajo 
estudio, siendo la T, la variable más dinámica 
y por ende la más limitante. En particular, se 
analizó la metodología propuesta por Jiang 
e Islam (2001), que permite la estimación 
directa de FE espacialmente distribuida, 
basándose en una interpretación contextual 
del espacio NDVI- T. Es por ello que en este 


trabajo se puso énfasis en analizar el efecto 
de la resolución espacial de las imágenes de 
T. y como éstas influyen en el espacio NDVI- 
T.. También se evaluaron estos efectos en el 
cálculo del parámetro 6 (Jiang e Islam, 2001), 
y por ende sobre FE, ya que esta última es 
directamente proporcional a (. 

El análisis estadístico de las imágenes de 
T. agregadas indicaría que la información 
térmica regional no cambia significativamente 
con la disminución de la resolución espacial, 
por lo que una imagen de T, a resolución 
espacial de 10 x 10 km puede representar 
la variación de temperatura de la zona en 
estudio. Sin embargo, la incertidumbre en 
las estimaciones hidrológicas derivadas debe 
analizarse cuidadosamente. 

Se estudiaron las características del dia- 
grama triangular, más específicamente del 
parámetro T.. de Jiang e Islam, y el borde 
inferior (representado por T. ). Las diferencias 
encontradas en el parámetro T',. son com- 
patibles con los errores comunes en el cálculo 
de T, (Kustas et al., 2003; Venturini et al., 2004). 
Las variaciones de T encontradas en este 
trabajo son comparables con las localizadas 
cuando se cambia de sensor, pero se mantiene 
la resolución de 1 km? o se trabaja con mejores 
resoluciones espaciales. Estos resultados 
indicarían que imágenes de T' con píxeles más 
grandes de 1 x 1 km no influirían en el borde 
superior del triángulo. 

Se destaca que durante los días 259/04 
y 266/04 se observan bordes fríos lineales 
no horizontales; sin embargo, los mismos se 
tomaron horizontales, tal como lo proponen 
los autores del método. Si bien se puede 
adoptar el borde frío como no horizontal, sin 
modificar los supuestos propuesto por Jiang e 
Islam, la determinación de la pendiente de la 
recta requeriría de más supuestos sobre dónde 
se asigna el valor de (, . Así, no se modificó 
el criterio de determinación del borde frío, a 
fin de no introducir incertidumbre al análisis, 
a pesar de lo observado. En dichos días se 
observaron las mayores deformaciones de los 
triángulos obtenidos y por consiguiente las 
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mayores diferencias en las imágenes de ( para 
las distintas resoluciones analizadas. 

Se obtuvieron los valores del parámetro 
con dos procedimientos diferentes. El prime- 
ro, denominado Ta, “onstante, que utiliza 
solamente el valor T,. calculado a partir de 
la imagen original con resolución 1 x 1 km. El 
segundo procedimiento consistió en repetir el 
cálculo de T, para de cada una de las imáge- 
nes agregadas. 

Al calcular Y con T,, constante, las 
diferencias encontradas en el valor medio 
regional son inferiores al 14% del valor medio 
de q . El método arrojaría resultados acepta- 
bles en términos de 4 para el área en estudio 
y todas las resoluciones espaciales de las 
imágenes de T, analizadas. Esto significa que 
el método da buenos resultados en la medida 
en que se rellenen datos faltantes con imágenes 
provenientes de MP, pero que los píxeles que 
determinan el T,, conserven la resolución de 
los sensores infrarrojos térmicos. 

Al estimar q con T', . variable, las diferencias 
encontradas en el valor medio de 4 llegan 
al 55% del valor medio de b, observándose 
las menores diferencias para imágenes con 
píxeles de 1x 1 y 2 x 2 km. Esto sugiere que 
el método arrojaría valores de q similares con 
imágenes cuya resolución espacial varíe entre 
1x1 y 4 x 4 km. Para resoluciones espaciales 
superiores, los resultados del método podrían 
diferir sensiblemente de los obtenidos con las 
imágenes a resolución original. 

En general se observa que el método es 
sensible al parámetro T', . En el caso de que 
dicho parámetro sea obtenido con imágenes 
de baja resolución espacial, por ejemplo con 
pixeles de 4 x 4 km, se observan diferencias 
importantes en relación con las imágenes de 1 
x 1 km. 

En general, el modelo desarrollado por 
Jiang e Islam sería aplicable con imágenes 
de baja resolución espacial siempre que se 
adopte el borde inferior del triángulo igual 
a la temperatura media del agua para una 
resolución de 1 x 1 km, independientemente de 
la resolución espacial de las imágenes de T, con 


las que se obtengan los espacios triangulares. 
Si no es posible calcular T', . con imágenes de 
T, con buena resolución espacial (por ejemplo, 
1 x 1 km), por no estar disponibles o porque 
la cobertura de nubes lo impide, las imágenes 
de T, derivadas de microondas pasivos solu- 
cionarían el problema parcialmente, teniendo 
en cuenta que los resultados de FE podrían 
tener diferencias de hasta el 50%. 
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Abstract 


GIROLIMETTO, D., VENTURINI, V. € RODRÍGUEZ, L. The effect of the spatial 
resolution of NDVET, spaces on the evaporative fraction calculation. Water Technology 
and Sciences, formerly Hydraulic engineering in Mexico (in Spanish). Vol. 1, No. 3, 
July-September, 2011, pp. 21-33. 


This paper analyzes the effect of the spatial resolution of the images used to generate the 
Normalized Difference Vegetation Index (NDVD-surface temperature (T) space. This 
is widely used to monitor environmental variables such as evapotranspiration (ET) and 
evaporative fraction (FE). The method proposed by Jiang and Islam (2001) was selected 
to study the effect of the spatial resolution on the calculation of the FE. For this, spatial 
resolutions were analyzed ranging from 1 x 1 km (common for thermal sensor images) to 10 
x 10 km (consistent with images from passive microwave sensors (MP)). The results suggest 
that the NDVIT, triangular space allows for obtaining reasonably accurate FE values with 
low spatial resolution, such as that of MP, provided that the cold edge can be determined 
using images with 1 x 1 km spatial resolution. Ifit is not possible to calculate the cold edge 
with images of 1 x 1 km resolution, the results would be uncertain and, therefore, should be 
carefully analyzed. 


Keywords: surface temperature, MODIS, evapotranspiration, NDVL 
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